X-ray and optical bursts and flares in YSOs: results from a 5-day 
XMU-Newton monitoring campaign of L1551 
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Abstract. We present the results of a five-day monitoring campaign with XMM-Newton of six X-ray bright young stellar 
objects (YSOs) in the star-forming complex LI 551 in Taurus. All stars present significant variability on the five-day time scale. 
Modulation of the light curve on time scales comparable with the star's rotational period appeared to be present in the case 
of one weak-lined T Tauri star. Significant spectral variations between the 2000 and the 2004 observations were detected in 
the (unresolved) classical T Tauri binary system XZ Tau: a hot plasma component which was present in the X-ray spectrum in 
2000 had significantly weakened in 2004. As XZ Tau N was undergoing a strong optical outburst in 2000, which had terminated 
since then, we speculate on the possible relationship between episodic, burst accretion, and X-ray heating. The transition object 
HL Tau underwent a strong flare with a complex temperature evolution, which is indicative of an event confined within a very 
large magnetic structure (few stellar radii), similar to the ones found in YSOs in the Orion Nebula Cluster. 
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1. Introduction 

Young stars are copious sources of X-ray emission; this, cou- 
pled with the fact that X-rays are not strongly affected by dust, 
makes X-ray observations a very useful tool in the study of 
the early phases of stellar evolution. In particular, X-ray obser- 
vations of young stars provide insights into high energy pro- 
cesses, such as coronal emission, and trace the stars' magnetic 
activity. While Class objects have not yet been conclusively 
detected in the X-ray, X-ray emission from Class I objects is 
relatively common; Class II, or classical T Tauri stars (CTTS) 
and Class III, or weak-lined T Tauri stars (WTTS), are al- 
most always intense X-ray sources (e.g.lPreibisch & Feig elsonl 
2005). 

While X-ray emission is interpreted in WTTS to be coronal 
in origin (a scaled-up version of the solar activity with similar 
energy production and emission mechanism), the phenomenol- 
ogy for CTTS and Class I objects is more complex than a sim- 
ple manifestation of enhanced solar-type magnetic phenom- 
ena, which suggests that accretion and the presence of rela- 
tively massive accretion disks may also play a role. Evidence 
of the role of accretion is provided by TW Hya, one of the 
nearest CTTS and thus one for which high-resolution X-ray 



Send offprint requests to: G. Giardino, 
Giovanna.Giardino@rssd.esa.int 



spectra are currently obtainable. Its inferred plasma temper- 
ature distribution, density, and peculiar chemical abundances 
are consistent with a model in which th e bulk of the X-ray 
emission is generated via mass accret ion jKastner et all l2002: 
IStelzer & SchmitA I2004I and see also lDrakell2005l) . The pres- 
ence of an accretion funnel shock at the site of the X-ray and 
UV emission is also invoked to explain the low flux in the for- 
bidden line in th e O vn triplet of the X MM-Newton RGS spec- 
trum of BP Tau JSchmittet all f2005). Nevertheless, the pres- 
ence of hot plasma, hotter than possibly generated by the accre- 
tion shock, shows that different X-ray emission processes must 
be going on at the same time in accreting YSOs. Statistically, 
high-accretion objects have lower Lx t han non-accreting ones 
(e.g. in the Orion Neb ula Cluster, ONC. lFlaccomio et allT2003l 
Preibisch etaill2005h . 

Temporal variability is a useful tool for distinguishing be- 
tween the underlying X-ray emission processes in different 
young stellar types. The observation, for example, of rotation- 
ally induced modulation of X-ray emission probes the spa- 
tial distribution of the emitting plasma and the time scales 
of its temporal evolution. Until recently, there were only a 
handful of reports of rotational modulation of X-ray emission 
from young main sequence stars: VXR4 5, a young fast rotato r 
star (P = 0.22 d) member of IC 2391 (Mari no et all 120031) : 
AB Dor, another fast rotator ZAMS K0 dwarf (P = 0.51 d) 



(Hussain et al., 2005); and EK Dra (a young sun analogue) 
with a period P = 2.7 d, for w hich indications of X-ray ro- 
tational modulation are given hv lGiidel et al l 99.4 . More re- 
cently, rotational modulation in the X-ray ligh t curve of sev- 
eral pre-main sequence stars (PMS) (Flacco mio et all Eo05) 
has been reported, based on the two-wee k Chandra monitorin g 
campaign of the ONC known as COUP JGetman et alll2005t) . 

In a previous 50 ks XMM-Newton observation of the star- 
forming complex LI 551 performed in 2000, Favata et al. 
(2003; hereafter FGM03) report interesting and unexpected 
variability in the CTTS XZ Tau, apparently neither due to 
'classical' flaring nor to simple rotational modulation. To fully 
understand the observed variability (whose timescales were 
clearly under-sampled by the 2000 observation) and the un- 
derlying processes, we performed a monitoring campaign of 
L1551 with XMM-Newton, composed of 11 exposures of 
roughly 9 ks each regularly spaced over 5 days, for a total 
integration time of about 100 ks. Although much more lim- 
ited in temporal coverage and total sensitivity than COUP, the 
L1551 campaign discussed here probed similar times scales, 
and the larger collecting area of XMM-Newton allows spec- 
tral variations to also be monitored on shorter time scales. 
Unfortunately the observations were affected (as discussed 
in detail later) by rather unfavorable background conditions, 
which made it impossible to fulfill all the original goals of the 
project. Nevertheless, we exploited the data to extract the max- 
imum amount of information possible on the variability of the 
seven bright X-ray sources in the field. 

The paper is organized as follows: the observations and the 
data reduction procedure are presented in Sect. |2] results for 
the individual sources are presented in Sect. [3] while a global 
analysis of their variability is presented in Sect. |4] The results 
are discussed in Sect.|5j while the conclusions are summarized 
in Sect.© 

2. Observations 

The XMM-Newton observation discussed in this paper con- 
sists of 11 exposures of roughly 9 ks each of the L1551 star- 
forming cloud. All observations were pointed in the same di- 
rection with the same roll angle, with bore-sight coordinates 
of RA 04:31:39, Dec 18:10:00. The exposures were regularly 
spaced over 5 days, starting March 4, 2004 at 16:02 UTC. All 
three EPIC cameras were active during the observations, in 
full-frame mode with the medium filters. The XMM-Newton 
EPIC cameras have a field of view of approximately 30 ar- 
cmin, with a PSF (in the center of the field of view) with a 
resolution of approximately 5 arcsec, and a useful bandpass of 
0.3 to 10 keV. The XMM-Newton mission is described e.g. in 
Ijansen et ai1 l l200ll) . 

Table \l\ summarizes the start time, time length, and live- 
time, in each camera, for each of the 11 exposure. For MOS1 
and MOS2, the exposure's typical length (elapsed time) is ap- 
proximately 8 ks, except for obs. 0501 and 0901, lasting 10.5 
and 13.7 ks, respectively. The PN exposures are usually shorter, 
lasting typically 6-7 ks, with a number of exceptions: obs. 0701 
and 0801 lasted only 2.5 and 3.4 ks, respectively, and obs. 0501 
and 0901 lasted 10.2 and 14.0 ks, respectively. PN is alive for 



a significantly smaller traction of the exposures duration than 
MOS1 and MOS2. Over the 1 1 exposures, the total observation 
time for PN is 82 ks, while for MOS1/MOS2, it is 107 ks, with 
corresponding livetimes of 68 ks and 105 ks. 

Figure ^ shows the light curves of the background (total 
counts from each camera at energies above 8 keV) in the three 
instruments. Many of the 1 1 exposures were significantly con- 
taminated by proton flares, in particular obs. 0201, 0501 0901, 
1001, and 1201 display a number of intense short duration 
flares and long duration (time scales comparable to the expo- 
sure times) episodes of high background. PN data are generally 
more contaminated than MOS1 and MOS2. The background 
level appears to modulated on time scales similar to the 2 day 
duration of XMM-Newton orbit. The background is high, as 
shown by comparison with the 2000 observation, for which 
only time intervals with less than 15 background counts per 
30 sec bin were retained. Applying such a strategy to this set 
of observations would result in the rejection of the whole 5-day 
campaign. The cleanest segments here have background levels 
that are at least twice as high, and the most contaminated ones 
reach background count rates up to 30 times higher. 

To recover the maximum amount of information possible, 
in spite of the high background, data were processed in two 
different ways, one aimed at determining time-resolved spec- 
tral parameters for the X-ray bright sources, the other aimed at 
recovering information for the weaker sources. The latter in- 
volved stacking, or merging, the 1 1 individual exposures into 
one single photon list. In both cases, only photons with energy 
ranging between 0.3 and 8 keV were retained. 

For the spectral analysis of the brighter sources, only time 
intervals affected by the strongest background spikes were re- 
moved, ensuring that only a small fraction of the three detec- 
tors' livetime was discarded from each exposure (never more 
than 10%). This implies that the background level in individ- 
ual exposures will be very different; while not fully satisfac- 
tory, this is the only strategy that allows us to maintain the time 
coverage of the original observation. The effectiveness of back- 
ground subtraction has been tested by verifying the consistency 
of the derived spectral parameters for different levels of back- 
ground filtering. 

The data were processed with the standard SAS V6.0.0 
pipeline. Source and background photons were extracted from 
the filtered event files using a set of scripts developed at 
Palermo Observatory. Source and background regions were de- 
fined interactively in each exposure separately using the ds9 
display software. The background was extracted from regions 
on the same CCD chip and at a similar off-axis angle to the 
source region. Response matrices ("rmf and arf files") appro- 
priate to the position and size of the source extraction regions 
were computed. 

The spectral analysis was performed using the xspec pack- 
age VI 1.2, after rebinning the source spectra to a minimum of 
20 source counts per (variable width) spectral bin. The spec- 
tral fits were all carried out in the energy range 0.3-7.5 keV, 
unless otherwise stated. All light curves and spectra discussed 
in this paper were background subtracted. Since V827 Tau and 




time [ s J 



Fig. 1. The time evolution of the background (total counts from 
the camera at energies E > 8 keV per 30 s bin) for the entire 
observation, separately plotted for each of the three detectors. 



V1075 Tau were unavailable in the two MOS cameras 1 and the 
MOS1 and MOS2 data in obs. 0901 were corrupted and un- 
usable, all light curves presented in the following section are 
from PN data. For the spectral analysis, we used PN data for 
all the sources except for V710 Tau, for which we used the re- 
sults from simultaneous fits to PN, MOS1, and MOS2 spectral 
data (for all exposures apart 0901) due to the source's weak- 
ness. For the 5 sources for which data from all three cameras 
are available, joint fits to PN, MOS1, and MOS2 spectral data 
yielded spectral parameters compatible with the ones derived 
from PN data alone. 

To properly compare the present observation with the 2000 
one, we also reprocessed the 2000 data with the new SAS 
pipeline V6.0.0. This resulted in a number of somewhat dis- 
turbing differences; in particular, for two sources (V827 Tau 
and V710 Tau, discussed in detail in Sect.[3j, for which appar- 
ently valid spectra resulted from the V5.3.1 pipeline used by 
FGM03, no 'good' photons are left in the results of the V6.0.0 
pipeline. Indeed, the fluxes derived by FGM03 for V827 Tau 
and V710 Tau are very different from the ones derived here 
for the 2004 observation, a discrepancy evidently caused by 
the problems in the V5.3.1 processing. For the other sources, 
the flux and spectral parameters derived in FGM03 and in the 
present work are generally in good agreement, although some 
important differences will be discussed later. 



1 V827 Tau fell out of MOS1 and was close to the camera edge in 
MOS2, V1075 fell out of MOS1 and, because of its weakness, it is 
hardly visible in some MOS2 exposures. 



z.7. me mergea data 

Most of the known X-ray sources in LI 551 are too weak to 
be visible in any of the eleven individual exposures. To try to 
detect fainter sources in the field, we stacked the individual 
exposures for each EPIC camera using the task merge in the 
SAS V6.0.0 pipeline. Prior to stacking, each exposure was fil- 
tered to reduce the background. To recover faint sources, the 
background must be reduced significantly below the high val- 
ues present in the data; the 1 1 exposures were filtered (before 
merging) at a threshold level of 100 background counts per 30 
sec bin for PN and 60 counts per 30 sec bin for MOS1 and 
MOS2. Table|2]summarizes the total (merged) live time before 
and after filtering. The PN camera was the one most affected 
by proton flares, and the restrictive filtering applied discards 
nearly half of the data. Thus, the remaining 'clean' (low back- 
ground) data in this observation is only 2/3 for PN and just 
over 4/3 for MOS1 and MOS2 of the clean time in the ob- 
servation analyzed by FGM03 (even though the total exposure 
is twice as long). The observation is therefore not deeper than 
the one of FGM03. We used the merged observation to derive 
integrated spectra of 5 of the 7 bright sources (the ones that 
did not undergo a flare during the observation) and compared 
the resulting spectral parameters with the average value of the 
spectral parameters derived from the individual exposures. The 
good agreement between the two is an important verification 
of the validity of the spectra derived from the single exposures, 
despite the high background levels. An X-ray image produced 
from the merged data is shown in Fig.|2] together with a /f-band 
2MASS image. 

Table 2. Filtering thresholds (background counts per 30 sec 
bin) and livetime before and after filtering for the merged data 
for the three detectors. See text for more details. 



Merged observations 


Filtering [Counts] 


Livetime [s] 


Instrument 


E > 8 keV 


Before After 


PN 


100 


68482 35053 


MOS1 


60 


89521 73483 


MOS2 


60 


89660 74994 



3. Analysis of the individual sources 

3.1. V826 Tau 

V826 Tau is a K7 SB2 WTTS llMundt et allll983l). The sepa- 
ration between the two comp onents is 0.06 A U Jjensen et all 
1994) and the period is 3 9 d JMathieu||l994 . A photometric 
period of 3.7 d is known (Bo uvier et allll995ri . 

During the 2000 XMM-Newton observations (FGM03), the 
count rate of V826 Tau increased by ^ 50% over the 50 
ksec. Also during this 5-day monitoring campaign, V826 Tau 
showed significant variability on different time scales (Fig.fJJ 
with a factor of ^ 2 amplitude. 

The spectral parameters for each exposure are summarized 
in Tabled All PN spectra could be fit with absorbed 2T plasma 



lable 1. Start time, exposure duration and hvetime tor the 1 1 individual exposures comprising the XMM-Newton campaign on 
L1551. 





Start time [yyyy-mm-dd hh-mm-ss UTC] 


Elapsed time [ks] 




Livetime [ks] 






PN 


MOS1 


MOS2 


PN 


MOS1 


MOS2 


r IN 


MOS1 


MOS2 


UZU 1 


2004-03-04 


16:25:28 


i f f\i nn 

16:03:09 


i f no no 

16:03:08 


7.0 


8.7 


8.7 




8.6 


8.6 


UjU 1 


2004-03-05 


03:40:49 


no 1 o o o 

03:18:32 


no 1 o o n 

03:18:29 


7.9 


9.6 


9.6 


/ . l 


9.4 


9.4 


U4U1 


~\ l \t \ \ / \ O IXC 

2004-03-05 


17:12:35 


16:50:24 


16:50:15 


7.0 


8.7 


8.7 


D.J 


8.5 


8.5 


UjUI 


o nn a no r\/~ 

2004-03-06 


01:03: 13 


nn /in n 

00:40:57 


nn /in r - o 

00:40:52 


10.2 


1 1.2 


1 1.2 


£ o 
O.y 


10.9 


10.9 


UoUl 


o r\r\ a no r\/~ 

2004-03-06 


16:17:29 


15:55:10 


15:55:08 


7.0 


8.7 


8.7 




8.6 


8.6 


U /Ui 


2004-03-07 


no ^ >i *o i 

03:44:21 


01:56:54 


01:56:54 


2.5 


9.3 


9.3 


z.z 


9.0 


9.0 


UoU 1 


'innd no m 

2004-03-07 


15:53:17 


14:46:39 


14:46:39 


3.4 


7.8 


7.7 




7.4 


7.4 




7004-03-08 


01:44:52 


01:22:32 


01 -22-32 


14.0 


15.6 


15.7 


1 n ^ 
1U. J 


15.2 


15.4 


1001 


2004-03-08 


16:19:29 


15:57:11 


15:57:09 


7.0 


8.7 


8.7 


6.3 


8.6 


8.6 


1101 


2004-03-09 


01:59:08 


01:36:55 


01:36:49 


8.7 


10.4 


10.4 


7.7 


10.1 


10.1 


1201 


2004-03-09 


14:53:25 


14:31:13 


14:31:05 


7.0 


8.7 


8.7 


6.1 


8.5 


8.5 





\: ■ . 



Fig. 2. The left panel shows a mosaic of 2MASS K-band images of the region observed by XMM-Newton, while the right panel 
shows an X-ray image produced with the merged low-background data from all three EPIC detectors. The two images are on the 
same scale and have been registered (field size ~ 30' x 30' centered at 04:31:44 +18:10:29). The sources discussed in the paper 
have been labeled. 



models with metallicity frozen at 0. 17 Z (the best-fit value ob- 
tained by FGM03). The count rate variability is not linked to 
large variations in the spectral parameters. The average values 
of the spectral parameters and of the X-ray luminosity are com- 
patible with the values determined by FGM03. 



The best-fit spectral parameters are reported in Table 
Excluding the flaring intervals, no significant spectral varia- 
tions were present. The typical values (N(H) ~ 0.1 x 10 22 cm -2 , 
kTi ~ 0.3 keV, kT 2 ~ 1 .2 keV, and Z = 0.14 Z ) are very simi- 
lar to the values derived in FGM03. 



3.2. V827 Tau 

V827 Tau is a K7 X-ray bright W TTS with a rotational period 
of 3.75 d jBouvieret all 1 19951) . In the light curve (Fig. 14} a 
large flare lasting over one day is present with a count rate in- 
crease of a factor of 10. The emission before the flare varies 
significantly, by a factor of 3 over ^ 3 d. 



As in Tabled the source intrinsic flux prior to the flare de- 
creases slowly by a factor of 3. The minimum X-ray luminosity 
is 4 times higher than the value reported by FGM03; however 
(as discussed in Sect.[2}, reprocessing with the V6.0.0 software 
has shown that the data for this source were not properly re- 
duced by the previous version of the SAS pipeline (probably 
because the source was very close to the edge of one of the 
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Fig. 3. Right panel - From top to bottom: light curve of V826 
Tau over the 5 days of monitoring (bin-time 3600 s), best-fit 
values for N(H) kT u kT 2 and for the ratio EMi/EM 2 . Left 
panel - light curve and best-fit values of spectral parameters for 
V826 Tau from the XMM-Newton observation of Sept. 2000 
(from FGM03). All data are from the PN camera. 



Fig. 4. From top to bottom: light curve of V827 Tau over the 
5 days of monitoring (bin-time 3600 s) and best-fit values of 
N(H) kT u kT 2 , and EMi/EM 2 ratio (PN data), for the nine ex- 
posures not affected by the star's flare. During the star's flare 
the parameter values of the absorbed 3T-plasma fit are given. 
The PN data from XMM-Newton observation of 2000 are not 
given for this source because they were found unreliable after 
reprocessing with more recent SAS pipeline (see text for more 
details). 



chips); thus the X-ray luminosity value of V827 Tau given in 
FGM03 is most likely incorrect. 



3.2.1 . Flare analysis 

To determine the flaring emission spectral parameters we sub- 
divided the rise phase of the flare (Obs. 1001) into two seg- 
ments (one for the rise phase and one for the peak). In order 
to separate the quiescent contribution from the flaring emis- 
sion, we modeled the spectra with an absorbed 3T-plasma. The 
parameters of the first two components (representing the qui- 
escent emission) were frozen to the values derived by simul- 
taneously fitting four spectra during the quiescent phase (Obs. 
0201, 0401, 0601, 0701, 0801). The values for the third com- 
ponent (the flaring emission) were then fitted. The results for 
this component are summarized in Table[3] 

To derive the flare's ph ysical parameters w e used the ap- 
proach initially discussed bv lReale et alJ (H"997) and since then 
applied to a variety of stellar flares. The calibration of the 
method for the XMM-Newton detectors, and a detai led expla- 
nation of the physics behind it can be found in iReale et all 
ll2004l) . to which the reader is referred. This approach allowed 
us to account properly for the presence of sustained heating 
during the flare decay, using the slope £ of the flare decay in 



the log T vs. log yEM diagram. The semi-length of the flaring 
loop in this formulation is given by 



tlc 



aF(0 



0.35 < £< 1.6 



(1) 



where a — 3.7 x 10 _4 cm _1 s _1 K 1/2 , tlc is the 1/e folding time 
of the light curve decay, and T m . dx is the peak temperature of 
the plasma in the flaring loop. Then, F(0 and the relationship 
between T m . dx and the best-fit peak temperature r b s are both 
functions that need to be separately determined for each X-ray 
detector, depending on its spectral response. For the EPIC PN, 



F(0 = cj(( - Z a ) + q a 



(2) 



where c a = 0.51 ± 0.03, £, = 0.35 ± 0.01, and q a = 1.36 ±0.18. 
The range of validity corresponds to an impulsively heated flare 
(£ — 1 .66) and to very slow decays (strong sustained heating) 
corresponding to the locus of statics loops (£ < 0.35). The max- 
imum temperature r max is derived from r o b s as 



o.i3or 



(3) 



The log T vs. log yEM diagram for the V827 Tau flare is 
shown in Fig. [5] Applying the above formalism to the flare, 
we derive tlc = 18.5 ks and a slope £ = 0.73 + 0.10. Applying 
Eqs.|2]andn] the resulting loop semi-length is L = 2.5-3.5,7? . 
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Fig. 5. The evolution of the V827 Tau flare in the log T vs. 
log y/EM plane. 



Given its bolometric luminosity and temperature (Lboi — 1-1L Q , 
T e s = 4060 K, as estimated by Bricen o et alll2002h . the radius 
of V827 Tau is R =s 2 R Q , so the loop has a size comparable to 
the star itself. Assuming for the loop a radius r — Q.1L, typical 
of solar events, from the emission measure at the flare maxi- 
mum one derives an electron density n e 1.3 x 10" crrT 3 , 
which is typi cal of the va lues found in the intense flares of 
YSOs (Favat a et all l200^) . The corresponding equipartition 
magnetic field strength is B ^ 340 G. 

3.3. V1 075 Tau 

VI 075 Tau is a K7 bin ary WTTS with a 2.43 d rotational period 
llBouvieretal.[ll995l) . In the 2000 observation the VI 075 Tau 
count rate decreased by a factor of ^ 2 in about 30 ks. The 
FGM03 spectral analysis showed that, while the temperature 
did not change significantly, the absorption varied from 0.19 + 
0.03 x 10 22 crrT 2 , when the source was more intense, to 0.08 ± 
0.03 x 10 22 cirT 2 at the lower flux level. 

The light curve of VI 075 Tau is shown in Figure [6] 
Variability over a range of time scales is present, including a 
flare in the second exposure and an apparent modulation on a 
time scale of a couple of days. 

The spectral parameters are reported in Table [8] Given the 
moderate statistics, the data could be satisfactorily fit with a 
single temperature plasma. No changes in N(H) are visible, 
while kT varies between 0.75 and 1.2 keV. Figure shows 
a scatter plot between the spectral temperature and the count 
rate. The two are clearly corre lated (P = 0.99993 from a 
Wilcoxon test, IWilcoxoii fl945). showing that the variability 
is more likely intrinsic, rather than due to rotational modula- 



Fig. 6. Right panel - From top to bottom: light curve of V1075 
Tau over the 5 days of monitoring (bin-time 3600 s) and best-fit 
values of N(H) and kT. Left panel - light curve and best-fit val- 
ues of spectral parameters for VI 075 from the XMM-Newton 
observation of Sept. 2000 (All are PN data). An absorbed 2T- 
plasma model was used for the spectral fits of the 2000 data. 
All data are from PN. 



tion (which would probably result in a 'gray' modulation, i.e. a 
spectral-independent change in the X-ray flux). 

The spectral parameters for the merged spectrum (N(H)= 
0.10 ± 0.01 x 10 22 cm" 2 , kT] = 0.37 ± 0.02, lcT 2 = 1.00 ± 0.02 
and Z = 0.19 ± 0.03) are similar to the FGM03 values, as is the 
X-ray luminosity (L x = 2.0 x 10 30 erg s -1 ). 

3.4. V710 TauA&B 

V710 Tau is a CTTS+WTTS binary JCarkner et allll996l) with 
spectral types Ml and M 3 and masses 0.68 Mq and 0.48 Mq 
(Jensen & Akeson, 2003). With a projected separation of 3.2 
arcsec, the system is not resolved by XMM-Newton. Its light 
curve is shown in Fig. [8] The system shows significant X-ray 
variability with the source counts varying by more than a factor 
of two. FGM03 found no evidence of variability from V710 
Tau during the 50 ks 2000 XMM-Newton observations. 

Due to the faintness of the source, PN and MOS data were 
fitted simultaneously with a single temperature model; Table [9] 
reports the spectral parameters. The merged PN spectrum of 
V710 Tau has sufficient statistics for a 2T fit. The spectral pa- 
rameters (iV(H)= 0.19 ± 0.06 x 10 22 cnr 2 , kT x = 0.30 ± 0.053 
kT 2 = 0.93 + 0.06, Z = 0.2 frozen) showed somewhat lower 
temperatures than in the 2000 observation (kTi = 0.63 + 0.05, 
kT 2 = 1.24 ± 0.10, FGM03). However, the SAS V6.0.0 re- 
processing of the FGM03 data flags most of the source pho- 



lable 3. Best-fit spectral parameters or the Haring component in V827 lau during Hare episode (Ubs. 1001 has been subdivided 
into 2 segments). Z = 0.3 Z (frozen); units are EM53 = 10 53 cm 4 . The parameters of the first two components were frozen 
to values derived by fitting four spectra during the quiescent phase (to the values N(H) = 0.09 x 10 22 cirT 2 , kT\ = 0.36 keV, 
EM 1 = 3.60 x 10 53 cm 3 , kT 2 = 1.16, EM 2 = 6.64 x 10 53 cm" 3 , and Z = 0.11 Z Q ) as described in the text. 
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Fig. 7. Count rate of VI 075 Tau versus plasma temperature for 
the 1 1 exposures. 



tons as invalid (due to a row of hot pixels). The reprocessed 
MOS2 data of V710 Tau are, on the other hand, 'clean', and the 
best-fit to the 2000 data gives N(H)= 0.30 + 0.08 x 10 22 cnr 2 , 
kTi = 0.38 ± 0.10, kT 2 = 1.12 ± 0.06 (P = 0.06), similar to the 
values of the current observation. The source X-ray luminosity 
from the merged data is 0.6 x 10 30 erg s _1 , a factor of 2 lower 
than the 2000 value (1.3 x 10 30 erg s" 1 ). 



3.5. XZ Tau 

XZ Tau is a bi nary CTTS with 0.3 arcsec separation 
llHaas et all Q])90), associated, together with HL Tau, with 
a co mplex set of bipo lar jets and Harbig Haro out- 
flows (Mun dtetall Il990l) . The spectral types are M2 and 
M3.5 for XZ Tau North and XZ Tau South, respectively 
( Hartigan & Kenvo nll2003l). A photometric period of 2.6 days 
has been derived bv iBouvier et al.l l[lj?95). They interpret this 
period as due to rotational modulation by a hot spot, 1500 K 



Fig. 8. From top to bottom: light curve of V710 Tau over the 
5 days of monitoring (bin-time 3600 s, PN data) and best-fit 
values of N(H) and kT (PN, MOS1 and MOS2 data). PN data 
from XMM-Newton observation of 2000 are not given for this 
source because they were found unreliable after reprocessing 
with more recent SAS pipeline (see text). 



hotter than the photosphere and covering 1.2% of the stellar 
surface. 

During the 2000 XMM-Newton observation (FGM03), the 
X-ray count rate increased by a factor of four in an approx- 
imately linear fashion over 50 ks. A time-resolved spectral 
analysis of the X-ray emission resulted in significant spectral 
changes, in particular a decrease in N(H) from 1 .06 x 10 22 cm -2 
to 0.26 x 10 22 cirT 2 . The temperatures increased from kT\ = 
0.14 keV and kT 2 = 2.29 keV to kT\ = 1.00 keV and kT 2 = 
4.98 keV (Table 5 in FGM03). The average spectrum was de- 
scribed well by a very low-metallicity plasma (Z = 0.007 Z Q ), 
with an X-ray luminosity Lx = 1.3 x 10 31 erg s" 1 . 

We reprocessed the data with SAS V6.0.0 and re-analyzed 
them, finding that the variation of N(H) reported by FGM03 
may be spurious. The re-analysis shows that the spectral data 



Fig. 9. Contours of x 2 as a function of the spectral parameters 
Af(H) and Z for the fit to the 2000 spectrum of XZ Tau repro- 
cessed with the V6.0.0 software, during the first time interval 
(the 50 ks observation was subdivided in three intervals). The 
X 1 space presents two well-separated minima that are both ac- 
ceptable solutions. The cross indicates the solution implying no 
change in N(H). 



for the first interval can be equivalently fit with two different 
solutions, one with the high N(H) reported by FGM03 and an- 
other with a low N(H) compatible with the later intervals. The 
two solution spaces are separated in the x 2 space well (as vis- 
ible in Fig- 0» and both have an acceptable best-fit probability. 
However, 'Occam razor' arguments lead us to prefer the one 
in which N(H) does not vary with respect to the rest. In fact, 
well-separated minima in the x 2 space are also present for the 
integrated 2000 spectrum. 

The results of the spectral fits to the reprocessed 2000 data 
are given in Table @] (replacing Table 5 of FGM03). The only 
spectral parameter that varies significantly over the 50 ks is 
the emission measure of the hotter component, which increases 
by a factor of five, with no appreciable increase in tempera- 
ture. A very-similar phenomenon - a large "slow" increase in 
the X-ray flux without significant change in the plasma tem- 
pera ture - was obse rved in the outbursting YSO VI 674 Ori 
JGrosso et alll2005i) . The average spectral parameters for the 
spectrum of the entire 50 ks observation resulting from this 
re-analysis are N(U) = 0.24 ± 0.01, kT x = 0.83 ± 0.02 keV 
kT 2 = 4.09 ± 0.34 keV, and Z = 0.22 ± 0.08 for a corresponding 
intrinsic luminosity of 2.3 X 10 30 erg s _1 , a factor of 6 lower 
than reported by FGM03, due to the change in N(H). 

The 2004 XMM-Newton light curve of XZ Tau is shown in 
Fig-C0l The data points relative to Obs. 0701 and Obs. 0801 
have been omitted as they are strongly contaminated from the 
nearby HL Tau, which during this time is undergoing a large 
flare (see Sect. |6]>. 

Table 1 1 01 summarizes the results of the spectral analysis of 
the PN data for the individual exposures; these were satisfacto- 
rily fit by an absorbed IT plasma model. The spectral fits to the 
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Fig. 10. Right panel - From top to bottom: light curve of XZ 
Tau over the 5 days of monitoring (bin-time 3600 s) and best-fit 
values of N(H) and kT (Obs. 0701 and Obs. 0801 have been). 
Left panel - light curve and best-fit values of spectral parame- 
ters for XZ Tau from the XMM-Newton observation of Sept. 
2000; for the spectral fits of the 2000 data an absorbed 2T- 
plasma model was used. All data are from PN. 



data for observations 0401 and 1101 can be improved by using 
a two-temperature plasma model (with larger error bars on the 
best-fit parameters); to allow a uniform comparison of the time 
variability we use the results of the IT spectral fits. 

Figure [Tolshows the best-fit values of N(H) and kT. While 
variability of a factor of ^ 2 in the count rate is present, no sig- 
nificant spectral variations are detected. A spectral fit to the 
merged data for XZ Tau (with the exclusion of Obs. 0701, 
0801, and 0901, to prevent contamination from the HL Tau 
flare) yields N(B)= 0. 19 ± 0.03 x 10 22 cm" 2 , kT x = 0.34 ± 0.03, 
kT 2 = 1.22 + 0.07, and Z = 0.14 ± 0.04 (P = 0.02, AE = 
0.3 - 3.0 keV) for a corresponding luminosity of 1.3 x 10 30 erg 
s~'(a factor of 2 lower than in 2000). 

The plasma temperature during the 2004 observation is 
significantly lower than in 2000, with no evidence of the hot 
(T a 4 keV) component clearly present in the 2000 data. 
Visual inspection of the two spectra (top panel of Fig. which 
plots both spectra reprocessed with the standard SAS V6.0.0 
pipeline) shows that indeed the high-energy tail clearly visible 
in the 2000 spectrum has weakened significantly. For compari- 
son, the bottom panel of Fig.[^shows the 2000 and 2004 spec- 
tra of V826 Tau, for which no differences are visible, allowing 
us to exclude e.g. problems in the background filtering. 

XZ Tau is the only star in the sample showing substan- 
tial spectral changes from 2000 to 2004. In 2000 XZ Tau 
had a hotter X-ray spectrum than in either 2001 or in 2004: 



lable 4. Best-fit spectral parameters tor XZ lau during three consecutive time intervals or the 2000 XMM-Newton observation. 
See the text for the details. 

Time interval N(H) kT x kT 2 £M, EM 2 Z X 2 p 

ks 10 22 cm' 2 keV keV 10" cm' 6 10" cm' 6 Z G 

0-20 0.22 ± 0.04 0.84 ± 0.03 4.31 ± 1.34 0.32 ± 0.59 0.49 ±0.10 0.40 ± 0.34 133 0.06 

20-40 0.22 ±0.021 0.82 ± 0.04 4.81 ±0.77 0.69 ± 0.83 1.97 ±0.57 0.20 ±0.12 0.79 0.95 

40-54 0.24 ±0.02 0.83 ±0.06 3.51 ±0.38 0.62 ± 0.90 2.46 ± 0.28 0.24 ±0.15 0.94 0.63 
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Fig. 11. Spectra of XZ Tau (top panel) and V826 Tau (bot- 
tom panel) from the merged PN data from the 2004 monitoring 
campaign (filled circles) and from the 2000 observations (open 
triangles). While the V826 Tau spectra are almost identical, the 
hard component present in 2000 in XZ Tau is no longer visible. 



the 2000 spectrum indicates kT\ = 0.83 ± 0.02 keV and 
kT 2 = 4.09 ± 0.34 keV, while the present (2004) data indicate 
kT\ = 0.3 - 0.7 keV and kT 2 = 1.0-1.5 keV, which are very 
similar to the values derived from the July 2001 Chandra ob- 
servation (FGM03), kTi = 0.65 + 0.03 and kT 2 = 1.56 ± 0.12. 
The light curve of the 2000 observation shows no obvious evi- 
dence of flaring (which would justify the higher temperature); 
slow variability is present, starting from a lower value than any 
of the present observations and rising to a value similar to the 
2004 one. 

Ground-based UVBR photometry of XZ Tau from 1962 to 

j 1 1 1 

1993 detected variations of > 2 mag (Herbst et al., 199 4t), inte- 
grated on both components of the binary system. ICoffev et alJ 
(2004]) monitored XZ Tau and its outflows from 1995 till 2001 
using HST, which resolves the binary. XZ Tau South (R 13.5) 
displays moderate variability (AR < 0.3 mag), while XZ Tau 
North (the suspected source of the outflow) displays dramatic 
variations: in Jan. 1995 its magnitude was R = 14.93, fading by 



about 1 mag until 1998 and thereafter brightening by 3 mag, 
so that by Feb. 2001 XZ Ta u North was a ctually the brighter 
star (also visible in Fig. 1 of lCoffev etail E004). This behav- 
ior suggests that XZ Tau North is an EXor. EXors, named after 
their prototype EX Lupi, are a loosely defined class of erup- 
tive CTTS that periodically undergo outbursts from the UV to 
the optical. Although increases by several magnitudes w ith ris e 
times of up to a few years have been recor ded fcferbiel ll989). 
the changes in these YSOs are not as extreme as in FU Ori stars 
(EXor spectra during outburst, for example, continue to resem- 
ble the ones of T Tauri stars). This phenomenon is thought to 
be due to major increases in the underlying disk accretion rate; 
however, the number of known EXors is relatively small and 
the class remains poorly defined. 

While the HST data presented bv lCoffev et"ail (12004 ex- 
tend only to Feb. 2001, a number of more recent HST observa- 
tions are present in the public archive, which can help to deter- 
mine whether the outburst that was ongoing at the time of the 
2000 observation has lasted into the 2001 Chandra and 2004 
XMM-Newton observations. We have located a number of un- 
published observations in the HST archive, which we extracted 
and inspected. While a full photometric analysis of these ob- 
servations is beyond the scope of the present work, even sim- 
ple visual inspection of the images allows assessment of the 
lumi nosity of XZ Tau N (the component found to strongly vary 
bv lCoffev et ail 12004 compared to the more stable XZ Tau S. 
We have chosen four HST observations (all taken in red filters) 
to be as close as possible in time to the X-ray observations. 
A zoom centered on the resolved XZ Tau binary is shown in 
Fig. El The first three observations are all from the WFPC2 
camera, in the F675W filter, while the last (2004) observation 
is taken with the ACS camera, with the FR656N filter. The 
XMM-Newton observation of FGM03 took place in Sept. 2000 
between the HST observation of Feb. 200 0, when the bright- 
ening of XZ Tau N was already going on (Coff ev et all 12004). 
and the HST observation of Feb. 2001, when the outburst was 
clearly visible. The outburst, however, did not appear to be long 
lasting: one year later (Feb. 2002) the outburst was already fin- 
ished, with XZ Tau N back to being fainter than the S com- 
ponent. In early 2004 it was, if anything, even fainter than the 
S component. The Chandra observation of July 2001 fell be- 
tween the early 2001 and early 2002 HST observations, and 
the present XMM-Newton campaign (4-9 Mar. 2004) is very 
close in time to the 2004 HST observation. 

While the X-ray and optical observations are not simulta- 
neous and the binary system is not resolved by XMM-Newton, 
the temporal proximity of the optical outburst of XZ Tau N 
and of the X-ray spectral hardening (as well as of the pecu- 
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Fig. 12. A mosaic of HST images of the XZ Tau binary taken 
between 2000 and 2004 (field size 4" x 3" centered on XZ 
Tau, at 04:31:40.07 +18:13:57.1). While during the outburst 
episode of 2001 XZ Tau N was the brightest component, 1 
year later it had faded back to being the fainter component, and 
by early 2004 (when the XMM-Newton observations discussed 
here have been taken) the S component is much stronger. The 
2000, 2001, and 2002 images are from the WFPC2 camera, in 
the F675W filter, while the 2004 is from the ACS camera in the 
FR656N filter. 



o.i 



o 

10 

>" 8 

E— 1 

* A 

2 




Fig. 13. Right panel - From top to bottom: light curve of HL 
Tau over the 5 days of monitoring (bin-time 3600 s) and best- 
fit values of N(H) and kT. Left panel - light curve and best- 
fit values of spectral parameters for HL Tau from the XMM- 
Newton observation of Sept. 2000. 



liar variability observed in 2000) is very suggestive of a con- 
nection between the two phenomena. S uch a connection ha s 
been well-documented for VI 647 Ori by Kastner et al. ( 2004), 
who reports a significant increase in spectral hardness near the 
peak of the optical outburst, as well as a surge of X-ray flux. 
Moreover, as already mentioned, the outbursting source VI 647 
Ori also displa ys the type of 'slow ' variability observed for XZ 
Tau in 2000 dGrosso et alll2005l) . 

3.6. HL Tau 

HL Tau, at about 24 arc sec from XZ Tau, is a K7 embedded 
young stellar object, often considered a prototype very young 
solar-mass star (0.5-0.7 Mq), with a circumstellar disk that re- 
sembles t he solar nebula at the early stages of planet formation 
dMen'shchikov et aTlll999h . Together with XZ Ta u, it is associ- 
ated with bipolar jets and Herbig-Haro outflows (Mundt et al., 
1990). 

Although initially classified as a CTTS, it is now thought 
to be in a transition phase from protostar (Class I object) to 
CTT S. An infall i ng en velope around the source was found 
bv lHavashi et al. (1993), who reported evidence of non-steady 
accretion. IClose et alJ (|1997) found that to reproduce the ob- 
served SED, the central source in HL Tau is required to be a 
very young (~ 10 5 yr) PMS surrounded by an active accretion 
disk and accreting at the rate of 5 x 10~ 6 Mq yr -1 . 

Ground-based photometry of HL Tau from 1973 to 1993 
has shown variations of 0.5 to 2 mag in UVBR dHerbst et all 



Il994l) . During the 2000 XMM-Newton 50ks observations 
HL Tau did not display signific ant variabi l ity, wh ile in the 2001 
80 ks Chandra observation (Bal ly et all 12003). it underwent 
a small short duration flare (Fig. 6 in FGM03). As shown in 
Figure^]in the present observation, HL Tau underwent a large 
(x 25 in count rate) flare which decays over about two days. 

The best-fit spectral parameters are reported in Table \H\ 
Given the low statistics of individual spectra, the metal abun- 
dance was frozen to Z = 0.6 Z , the value determined in 2000 
by FGM03. Outside of the flare the spectral parameters do not 
vary significantly and their average values (see below) are sim- 
ilar to the 2000 values (FGM03), as is the X-ray luminosity 
(1.5 x lO^ergs- 1 ). 

To derive the flaring emission spectral parameters, we ap- 
plied the same procedure as for V827 Tau. The quiescent spec- 
tral parameters are N(H)= 2.43 x 10 22 cm" 2 , kT = 3.10 keV, 
EM = 1 .40 x 10 53 cm 3 , and Z = 0.6 Z , and the resulting spec- 
tral parameters for the flaring component are listed in Table [5] 
The event shows a peculiar evolution, with a monotonically de- 
caying light curve associated with a highly irregular temper- 
ature evolution. The temperature has two well-defined peaks 
above 7 keV separated by a deep minimum at about 3 keV. 
The temperature evolution suggests that this is the combina- 
tion of two flares, probably physically related to each other but 
occurring in independent coronal structures. This kind of evo- 
lutio n has been predicted by modeling tw o independent flares 
hv lReale. Giidel. Peres. & AiidardN2004l) . 



lable 5. Best-ht spectral parameters or the flaring component in HL lau. Units are N22 = 10 cm , = 10 cm and 

F_i3 = 10~ 13 erg cm -2 s _1 . 
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In order to support this hypothesis, we modeled the event 
by combining two flares computed with detailed hydrodynamic 
modeling of plasma confined in a coronal loop. The light curve 
decay time suggests long flaring structures (e.g. ISerio et all 
Il99ll) . so that we considered each model flare to be identical to 
the one used to describ e one of the fla res observed during the 
COUP campaign (Fava ta et all 12005) in detail. We assumed 
that each flare occurs in a coronal loop with a constant cross- 
section and half-length L = 10 12 cm, symmetric around the 
loop apex. Both flares were triggered by injecting a heat pulse 
in the loop, which was initially at a temperature of 20 MK. 
This heat pulse is symmetrically deposited at the loop foot- 
points with a Gaussian spatial distribution of intensity 10 erg 
cm -3 s -1 and width 10 10 cm (1/100 of the loop half-length). 
After 20 ks the heat pulse was switched off completely. From 
the evolution of the plasma density and temperature along the 
loop co mputed with the Palermo-Harvard h ydrodynamic loop 
model ( Per es et al l ll982tlBetta etall 1997), we synthesized 
the corresponding EPIC spectra of the loop throughout the 
flare, deriving a light curve and the evolution of temperature. 

To model the HL Tau event we duplicated, the resulting 
light curve and temperature evolution with a time shift. The 
two flares are identical flares, except for a normalization fac- 
tor, which represents the loop cross-section and does not enter 
explicitly in the hydrodynamic modeling. The second flare has 
a normalization factor of 0.3 (i.e. a correspondingly smaller 
cross-section) and it starts 60 ks after the first. We summed 
the resulting two asynchronous sequences of flare spectra and 
obtained a single sequence of spectra, which we integrated to 
derive a single light curve and fit with single temperature EPIC 
model spectra. Figure ^] shows the resulting light curve and 
temperature evolution as compared to those obtained from the 
data. The model temperatures in the first flare are somewhat 
higher than the observed one, but the main flare characteristics, 
i.e. the monotonic light curve and the temperature dip, are re- 
produced by the double-flare model well - although the model 
is not unique. 

Since no constraint can be derived from the data, each of 
the two flares was modeled with no signific ant re s idual h eating 
present during the flare decay (e.g. iReale et all 119971) . Such 
modeling implies very large flaring structures, similar to the 
ones found in ONC YSOs bv lFavata etail J2005ll - where the 
data allowed investigation of the presence of sustained heating. 
Such large structures, with L 5R„, have only bee n found in 
YSOs, and were interpreted by iFavata et al.l J2005I) as linking 
the star to the accretion disk, i.e. as being the magnetic struc- 
tures supporting the magnetospheric accretion. In addition to 
the evidence from the ONC YSOs, HL Tau is the first Taurus 



YSO in which such large flaring structures have been detected. 
We cannot a priori exclude that shorter loops with sustained 
heating in the decay may also reproduce the features of this 
flare. However, the long delay of the model flares required to re- 
produce the distant temperature peaks suggests that very l arge 
structures must be involved in the flare (Real e et alll20Q4 . 

3.7. HD 285845 

Unlike the other stars in the present study, HD 285845 (an ac- 
tive binary system) is not a member of the star-forming re- 
gion, on the basis of its radial velocity and proper motion 
( Walter et al., 1988). The primary spectral type is G8, and 
Schneide r et alJ Jl998) report a separation of 73 mas and a 
magnitude difference of 1.19 mag. 

In the 2000 observation, HD 285845 showed significant 
variability, and similar behavior is present in the present data 
set (Fig.ll5>. Tablell2lsumm arizes the best-fit spectral parame- 
ters. Notwithstanding the significant count-rate variability, the 
spectral parameters show little variation in time. 

4. Temporal variability 

All the stars in the sample show significant variability, although 
with different characteristics. To quantify this variability we 
computed the normalized cumulative distributions of the am- 
plitude variability. These represent the fraction of time that a 
source spends in a state with the flux larger than a given value, 
expressed in terms of a given a normalization value, which can 
be the minimum count rate, the median count rate, etc. For the 
present sample we took as normalization value the count rate 
above which a source spends 90% of the time; this is less sen- 
sitive to noise fluctuations than the real minimum. 

Figure[H)]shows the distribution for the stars in our sample 
(except V710 Tau for which the count rate uncertainties are 
too large). Some sources (V826 Tau, HD285845 and XZ Tau) 
show mostly low-amplitude variability, in which less than 30% 
of the time is spent in a state 1.3—1.5 above the minimum. On 
the other hand, V827 Tau and HL Tau spend more than 60% of 
the time in such a state, while VI 075 Tau shows intermediate 
behavior. The differences between V827 Tau and HL Tau are 
not due to the large flares present in their light curves, as the 
flare points affect only the 10 bins relative to the highest count 
rate, while the low-variability tails of the distributions are also 
significantly different. 

Kolmogorov-Smirnov tests for the cumulative distributions 
of various pairs of sources also indicate the presence of the dif- 
ferent behaviors. While this would suggest differences in the 



1.00 



CD 

<5 



o 
O 



0.10 - 



0.01 



200 




150 



150 - 



100 - 




150 



Time [ks] 



Fig. 14. Top panel - The observed light curve for HL Tau (as- 
terisks), together with the light curve predicted by the flare 
modeling described in the text. Bottom panel - The observed 
temperature evolution of the HL Tau flare, together with the 
double-peaked temperature evolution predicted by the same 
flare modeling. The model is obtained by combining two mod- 
els of 2 x 10 12 cm-long flaring loops (dotted and dashed lines). 
The two models are identical, except for the cross-section of 
the flaring loops and for the start time (60 ks delay). 
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Fig. 15. Right panel - From top to bottom: light curve of 
HD285845 over the 5 days of monitoring (bin-time 3600 s), 
best-fit values of N(H) kT x , kT 2 and EM 1 /EM 2 ratio. Left panel 
- light curve and best-fit values of spectral parameters from the 
XMM-iVwton observation of Sept. 2000 (All are PN data). 
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X-ray emission processes, the two groups are heterogeneous, 
and both contain CTTS and WTTS, so that a simple interpreta- 
tion in terms of e.g. accreting vs. non-accreting sources is not 
possible. 



Fig. 16. Normalized cumulative count-rate distribution for all 
the star considered apart from V710 Tau. 



5. Discussion 

The X-ray monitoring observations of the LI 551 star-forming 
complex discussed here, in conjunction with the 2000 XMM- 
Newton observations (FG M03) and the 2001 Chandra obser- 
vations (Bally et al., 2003), has allowed us to probe the X-ray 
intensity and spectral variability from six YSOs (both CTTS 
and WTTS), on a variety of time scales, from hours to days to 
years. 

Most stars in the sample show X-ray variability at an am- 
plitude of a factor of ^ 2 on time scales of a few days. In the 
WTTS V1075 Tau, the source count rate varies on time scales 
similar to its optically derived rotational period, which would 
be compati ble with rotational mod ulation (as observed in YSO 
in the ONC. lFlaccomio et alll2005l) . but the significant correla- 
tion with spectral variability would appear to point to intrinsic 
variations. Plasma temperature variations of the order of 20% 
on similar time scales are also present, and in some cases they 
appear to be correlated with the intensity (e.g. for V1075 Tau). 

Again XZ Tau emerges as a peculiar source. While the 
short-term (day scale) spectral variability of XZ Tau is not 
linked to variations in the absorbing column density, as origi- 
nally speculated by FGM03 - reprocessing and reanalysing the 
original 2000 XMM-Newton data showing that the variation 
reported by FGM03 was probably spurious - this is the only 
source in the sample to show significant long-term changes in 
its spectrum. 

The spectrum of XZ Tau observed with XMM-Newton in 
2000 had a significant hard component (T =s 4 keV), which 
was not visible in the 2001 Chandra data and is not present 
in the 2004 XMM-Newton data and both are compatible with 
a hard component of at most 1.5 keV. The difference between 
the two spectral states is very well visible in Fig. ^2 The X-ray 
spectral hardening takes place close in time to the strong optical 
outburst of the N c omponent of the binary. While the data by 
ICoffev et alJ fc004 . who reported the outburst, only cover its 
beginning, our own inspection of new HST data indicates that 
the outburst was not long-lived, with the source returning to its 
state fainter than XZ Tau S by early 2002, one year after the 
early 2001 peak. 

Finally, while (Fig. llOt the count rate varied smoothly be- 
tween 0.04 and 0.22 cts/sec during the 2000 observation, during 
the 2004 campaign it never decreased below 0.12 cts/sec (with 
a peak value of 0.35 cts/sec), again showing the source to be 
in a different state. 

The binary system XZ Tau is not resolved in the X-ray, 
therefore one cannot exclude that XZ Tau S may have been 
the dominant X-ray source, at one or all the epochs at which it 
was observed. Nevertheless the observed phenomena are very 
reminiscent of the spectral hardening and variability observed 
in the outbursting source, VI 647 Ori iKastner et all [2004; 
iGrosso et alll2005l) . so that it seems reasonable to interpret the 
X-ray spectral changes and variability of XZ Tau as connected 
with the outburst of XZ Tau N. The potential physical mech- 
anism causing the significantly harder X-ray spectrum (with a 
somewhat lower emission level) during the optical outburst is 
difficult to assess. While the optical outburst is naturally inter- 
preted as an accretion outburst (with the optical luminosity due 



to increased luminosity of the accretion shock), why this should 
cause a hot plasma component to appear in the star is not clear. 
Indeed, accreting YSOs in the COUP sample are statistically 
less X -ray luminous than non-accreting ones llPreibisch et all 
2005) an d less prone to show flares from large magnetic struc- 
tures jFavataet alll2005l) . 

The temperature of the hot component observed in XZ Tau 
in 2000 is such that it cannot be due to simple accretion-driven 
emission, as the shock temperature in a low-mass T Tau star is 
too low . Its early M spectral type corresponds, at this age, to 
a typical mass M ^ 0.5 M and to a radius R st R Q , resulting 
ICalvet & GullbrineU 19981) in a peak shock temperature T - 
0.2 keV. 

Nevertheless, a heating of the plasma during an accretion 
outburst in XZ Tau N would be consistent with the observed 
high X-ray plasma temperatures in accre ting objects, like Class 
I sources. In a study of the p Oph region, Imani shi et all J2003) 
find that Class I sources have a higher (kT) (sometimes exceed- 
ing 5 keV) on average than Cl ass II-III sources. Th is is con- 
firmed by the study of p Oph bv lOzawa et all J2005I) that finds 
an evolutionary trend from Class I sources showing higher tem- 
perature and larger absorption to Class II and III sources show- 
ing lower temperatures and smaller extinction. In the present 
sample the only star with a consistently high plasma temper- 
ature (kT = 2.5-4 keV) is HL Tau, a PMS in the transition 
phase between Class I and Class II, which is still undergoing 
substantial accretion. 

Apart from XZ Tau, no significant variations in the activity 
level of our target YSOs on time scales of four years is ob- 
served. Whi le evidence of solar-like cyclical variati on is ac- 
cumulating JFavata et all 120041 iRobrade et all EooH) for low- 
and intermediate-activity stars, long-term variations (whether 
cyclical or not) in the activity levels of high-activity stars (in- 
cluding YSOs) do not seem to be present. 

Although the data do not allow a fully detailed model- 
ing, the double-temperature peaked flaring event observed in 
HL Tau most likely originates in very large magnetic struc- 
tures, similar to the ones observed in ONC YSOs. Such struc- 
tures, postulated by the magnetospheric accretion model, prob- 
ably link the stellar photosphere to the inner rim of the accre- 
tion disk. In the ONC sample only objects that are currently 
not accreti ng were shown to have very large flaring magnetic 
structures JFavata et alll2005l) . The only active accretor present 
in the sample showed flares confined to small structures, sim- 
ilar to normal coronal flares. HL Tau (with the caveat linked 
to the limited spectral information available) would provide a 
counter-example to the scheme observed in the ONC, being the 
first strongly accreting YSO with large flaring magnetic struc- 
tures. 



6. Conclusions 

The present monitoring campaign has allowed us to study the 
luminosity and spectral variability of the X-ray emission from 
a number of YSO over a range of time intervals, sampling in 
particular the variability over a 5-day span, as well as over 4 
years. 



While most YSOs show a remarkably constant level or ac- 
tivity over the 4-year span sampled, XZ Tau has shown signif- 
icant spectral variability apparently in conjunction with an op- 
tical outburst of the N component of the binary system, which 
took place in 2000. The role of accretion in the X-ray of YSOs 
is a matter of current debate, with a few spectral observations 
indicating differences in the density and UV environment of the 
plasma in a few accreting YSOs, and with a number of authors 
speculating on t he possible relationsh ip b etween accretion and 
coronal heating. Kastner et al. (2004) and lCrosso et'aTI ll2005l) 
report on spectral hardening and enhanced variability in VI 647 
Ori in conjunction with a well-documented optical/IR outburst, 
which provides evidence that strongly enhanced high-energy 
emission can occur as a consequence of high accretion rates. 
The significant spectral changes detected here supply similar 
evidence. 
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lable 6. Best-tit spectral parameters ot the FN data tor the eleven exposures or V826 lau. EM is the emission measure, f the 
null-hypothesis probability of the fit,and Fx Intr. the "unabsorbed" X-ray flux. Units are N22 = 10 22 cm~ 2 , EM53 = 10 53 cm -3 , 
and Fin = 10~ 13 erg cm 2 s _1 . The metallicity was frozen atZ = 0.17 Z©. The spectral fits were carried out in the energy range 
0.3-7.5 keV. 
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Table 7. Best-fit spectral parameters of the PN data for the eleven exposures of V827 Tau. EM is the emission measure, P the 
null-hypothesis probability of the fit, and F x Intr. the "unabsorbed" X-ray flux. 

Units are N22 = 10 22 crrT 2 , EM 53 = 10 53 crrT 3 , and f_i 3 = 10~ 13 erg cirT 2 s~'.The spectral fits were carried out in the energy range 

0.3-7.5 keV. 
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Table 8. Best-fit spectral parameters of the PN data for the eleven exposures of V1075 Tau. In this case Z was frozen at 0.15 Z Q , 
the value derived for this source from the XMM-Newton observation of L155 1 of 2000 (FGM03). EM is the emission measure, P 
the null-hypothesis probability of the fit, and F x Intr. the "unabsorbed" X-ray flux. Units are N22 = 10 22 crrT 2 , EM^ = 10 53 crrT 3 , 
and f _o = 10~ 13 erg crrT 2 s _1 .The spectral fits were carried out in the energy range 0.3-7.5 keV. 
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lable y. Spectral parameters derived from simultaneously htting FN, MOM, and MOS2 data tor the eleven exposures or V710 
Tau. In this case Z was frozen at 0.2 Z . EM is the emission measure, P the null-hypothesis probability of the fit, and Fx Intr. the 
"unabsorbed" X-ray flux. Units are N22 = 10 22 cirr 2 , £M 53 = 10 53 cm -3 , and F_i 3 = 10~ 13 erg cm -2 s _1 . The spectral fits were 
carried out in the energy range 0.3-7.5 keV. 
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Table 10. Best-fit spectral parameters of the PN data for the nine exposures of XZ Tau. In this case Z was frozen at 0.08 Z . 
EM is the emission measure, P the null-hypothesis probability of the fit, and Fx Intr. the "unabsorbed" X-ray flux. Units are 
N22 = 10 22 cm' 2 , EM53 = 10 53 crrT 3 , and F_i 3 = 10~ 13 erg crrT 2 s _1 . The spectral fits were carried out in the energy range 
0.3-7.5 keV. The quality of data was such that the parameter space of the model was degenerate in N(H) and Z, and changing the 
initial values could result in significantly different best-fit values for these two parameters. The value Z = 0.08 Z was determined 
by first performing a fit in which Z was allowed to vary, determining the average value over the 1 1 exposures and then repeating 
the spectral fits with Z frozen. 
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Table 11. Best-fit spectral parameters of the PN data for the eleven exposures of HL Tau. In this case Z was frozen at 0.6 Z , 
the value derived for this source from the XMM-Newton observation of L1551 of 2000 (FGM03). Spectral fits were limited to 
energy range 1 .0-7.5 keV. EM is the emission measure, P the null-hypothesis probability of the fit, and Fx Intr. the "unabsorbed" 
X-ray flux. Units are N22 = 10 22 crrT 2 , £M 53 = 10 53 crrT 3 , and F_ 13 = 10~ 13 erg crrT 2 s _1 . 
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2.73 ± 1.16 


2.36 ± 2.29 


0.5 


0.89 


3.8 


7.8 


0301 


2.19 ±0.58 


4± 1.78 


0.79 ± 0.56 


0.89 


0.51 


2.1 


3.7 


0401 


1.47 ± 0.8 


5.41 ± 3.33 


0.96 ± 0.76 


1.45 


0.14 


3.4 


4.4 


0501 


3.98 ± 1.94 


2.48 ± 1.46 


1.59 ±2.57 


0.64 


0.84 


2.2 


6.0 


0601 


3.04 ± 0.44 


1.75 ± 0.27 


2.66 ± 1.23 


0.7 


0.73 


2.7 


8.9 


0701 


3.14 ±0.19 


9.1 ± 1.83 


14.08 ± 1.03 


0.75 


0.9 


48.4 


81.9 


0801 


3.14 ±0.23 


4.56 ± 0.67 


12.75 ± 1.75 


1.03 


0.42 


33.2 


61.5 


0901 


3.41 ± 0.2 


8.53 ± 0.96 


8.46 ± 0.70 


0.93 


0.71 


28.0 


48.9 


1001 


2.87 ± 0.53 


2.41 ±0.53 


3.41 ± 1.73 


0.59 


0.95 


5.2 


12.7 


1101 


2.12 ±0.39 


4.11 ± 1.22 


1.34 ± 0.63 


1.49 


0.11 


3.7 


6.3 


1201 


2.08 ± 0.49 


3.4+ 1.16 


1.27 ±0.79 


1.65 


0.09 


3.1 


5.6 



lable 12. Best-tit spectral parameters or the bFlC-FN data tor the eleven exposures or HDzo5s45. EM is the emission measure, f 
the null-hypothesis probability of the fit, and Fx Intr. the "unabsorbed" X-ray flux. Units are N22 = 10 22 cirT 2 , £M 53 = 10 53 cm -3 , 
and F_i3 = 10~ 13 erg cm -2 s _I . The metallicity was frozen at Z = 0.14 Z©. The spectral fits were carried out in the energy range 
0.3-7.5 keV. 



Obs. 


N(H) 


kTi 


EM\ 




kT 2 


EM 2 


2 

X 


P 


fx 


F x Intr. 




N 22 


keV 


EM 53 




keV 


EM S3 






F-n 


F-n 


0201 


0.07 ± 0.02 


0.41 ± 0.06 


2.51 ± 1 


34 


0.97 ± 0.03 


4.94 ± 0.55 


1.07 


0.29 


15.16 


20.99 


0301 


0.11 ± 0.02 


0.34 ± 0.02 


4.79 ± 1 


39 


0.96 ± 0.03 


3.98 ± 0.39 


1.07 


0.27 


12.77 


21.19 


0401 


0.09 ± 0.01 


0.34 ± 0.02 


4.71 ± 1 


34 


1.00 ± 0.03 


5.19 ± 0.42 


1.04 


0.33 


16.13 


25.02 


0501 


0.11 ± 0.02 


0.36 ± 0.04 


4.29 ± 1 


84 


1.02 ± 0.04 


3.64 ± 0.54 


0.98 


0.57 


12.08 


19.77 


0601 


0.07 ± 0.01 


0.34 ± 0.03 


3.45 ± 1 


16 


0.98 ± 0.02 


5.57 ± 0.41 


0.99 


0.52 


16.87 


23.96 


0701 


0.07 ± 0.02 


0.35 ± 0.05 


4.04 ± 1 


96 


1.19 ± 0.05 


6.26 ± 0.61 


1.08 


0.3 


20.71 


28.09 


UOUl 


n HQ -t- n m 


V.JJ ± V.VJ 


^ HQ -i- 9 
J.uy ± Z 




1 .uz ± u.uo 




u.o / 


U. / D 


1Z. / J 


IXJXjy 


0901 


0.08 ± 0.02 


0.34 ± 0.03 


4.01 ± 1 


39 


1.01 ±0.03 


4.86 ± 0.47 


1.14 


0.03 


15.46 


22.71 


1001 


0.07 ± 0.01 


0.39 ± 0.03 


4.11 ± 1 


39 


1.19 ±0.04 


5.42 ± 0.44 


1.27 


0.01 


19.03 


26.05 


1101 


0.08 ± 0.01 


0.32 ± 0.02 


3.81 ± 1 


05 


1.05 ±0.03 


4.35 ± 0.33 


0.85 


0.93 


13.81 


20.58 


1201 


0.08 ± 0.01 


0.42 ± 0.06 


2.33 ± 1 


13 


0.97 ± 0.03 


5.33 ± 0.45 


0.99 


0.52 


15.56 


21.99 



